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Abstract 

We study instanton solutions and superpotentials for the large number of vacua 
of the plane-wave matrix model and a 2+1 dimensional Super Yang-Mills theory on 
Rx with sixteen supercharges. We get the superpotential in the weak coupling 
limit from the gauge theory description. We study the gravity description of these 
instantons. Perturbatively with respect to a background, they are Euclidean branes 
wrapping cycles in the dual gravity background. Moreover, the superpotential can 
be given by the energy of the electric charge system characterizing each vacuum. 
These charges are interpreted as the eigenvalues of matrices from a reduction for the 
1/8 BPS sector of the gauge theories. We also discuss qualitatively the emergence 
of the extra spatial dimensions appeared on the gravity side. 
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1 Introduction 



Many supersymmetric field theories have a large number of supersymmetric vacua. In 
the context of AdS/CFT correspondence P,[2],[3], when these field theories have string 
theory dual descriptions, these vacua are dual to supergravity backgrounds. These back- 
ground geometries are the saddle points where string perturbation theory can be ex- 
panded around. 

In this paper we study instantons interpolating between different supersymmetric 
geometries corresponding to different vacua of the gauge theories. We associate each BPS 
vacuum geometry with a value of a superpotetial, which is the vacuum expectation value 
of a particular operator on the gauge theory side. The superponetial can be calculated 
independently from gauge theory and gravity side. The superpotential from the gravity 
side reveals the quantum nature of the emergent geometry [4] (see also [5], [6], [7]) as a 
strong coupling effect of gauge theory. The emergence of extra spatial dimensions, extra 
relative to the boundary gauge theory, can be explained in terms of the emergence of the 
geometries from matrices in these gauge theories. 

The specific gauge theories we will study for these issues are the 0+1 dimensional 
plane-wave matrix model [8j and a 2+1 dimensional super- Yang-Mills theory (SYM) on 
R X 13 with sixteen supercharges. We chose these theories for study since the 
explicit formulations of these theories on gauge theory side and their supergravity dual 
backgrounds on the gravity side are explicitly available [10], [7]. The fluctuations around 
these backgrounds are dual to perturbations around the vacua of these gauge theories. 
These theories can be viewed as reductions of the 3+1 dimensional A/" = 4 SYM on 
and on the Hopf fiber of or S^/Zk. As a result, the string spectrum around these vacua 
share some similarities especially in the 5*0(6) sector [10]. In [TO], the near-BPS string 
spectrum was studied and shows a general strong-weak coupling problem. 

In order to better understand the vacuum structures of these theories, we study the 
effects of instantons interpolating between different vacua. On the gauge theory side, we 
analyze instanton solutions and instanton actions in the weak coupling regime. We find 
that it is convenient to define a superpotential and the difference of the superpotential 
between two vacua gives the instanton action. In the gravity side, we try to find the dual 
description of these instantons. We find in a large class of situations, they are described 
by Euclidean branes wrapping some cycles in the BPS geometries. The condition for 
instanton solutions in these cases is described by the condition of embeddings of these 
branes. We compute their actions and these give the strong coupling results of the 
instanton actions. These instantons interpolate between initial and final vacua which are 
very close to each other. 

To understand the instanton action corresponding to initial and final vacua that are 
not close to each other, we need to understand the instanton beyond the Euclidean brane 
approximation on a background. We need to know the non-perturbative superpotential of 
each vacuum. Each vacuum in the strong coupling regime is described by an electrostatic 
configuration with charges distributing on the disks under an external potential. The 
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charges of the electrostatic system can be interpreted as the eigenvalues of the matrices 
from the gauge theory side. We found that the superpotential for each vacuum is given 
by the energy of the eigenvalue system. We found consistency and agreement of the 
superpotentials between gauge theory and gravity side, except that the gravity result 
is the strong coupling result while the result from classical gauge theory action is in 
the weak coupling limit. As the couplings increase the actions of these theories receive 
quantum corrections. The gravity answer shows quantum repulsions of the eigenvalues, 
which was absent in the weak coupling limit, a feature similar to the Dijkgraaf-Vafa 
matrix model 

The non-perturbative superpotential is intimately related to the emergence of the 
geometries. We then study the emergence of all the extra spatial dimensions, relative to 
the boundary gauge theories. One of our methods is to embed the sector of these vacuum 
geometries into a larger sector of 1/8 BPS geometries associated with three harmonic 
oscillators, similar to [4j. We also try to explain the emergence of the electrostatic 
system and the forces between eigenvalues from gauge theory. On the gravity side the 
forces are described by the differential equations governing the gravity solutions. We 
speculate on other possible approaches. 

The organization of this paper is as follows: In section 2, we study instanton solutions 
and superpotentials of these vacua on gauge theory side in the weak coupling regime. The 
instanton actions are given by the differences of the superpotentials between initial and 
final vacua. In section 3, we study descriptions of these instantons and superpotentials 
on gravity side. These include Euclidean brane analysis in the cases when initial and final 
vacua are close to each other, and then an electrostatic description of the superpotential 
for general vacua. These are results in strong coupling regime of the gauge theories. In 
section 4, we discuss the emergence of the extra coordinates and the electrostatic system 
on gravity side, which encodes the information of the non-perturbative superpotential in 
the strong coupling regime. Finally in section 5, we conclude and discuss some related 
issues. 



2 Instantons on the gauge theory side 

The U{N) plane wave matrix model has a gauge field Aq, three S0{3) scalars Xj, i = 
1,2,3 and six S0{6) scalars X^, a = 4, 5,..., 9 and their fermionic partners. They are 
all X X X Hermitian matrices. The action is written in for example [8] . [14] . [T3] ■ [15] . [T2] . 
Here we use the action such that it depends only on a dimensionless coupling constant 
Qy^Q/m^, where m is originally the mass of the 5*0(6) scalars, and can be set as 1 for 
convenience. We have rescaled all fields so that the masses of 5*0(6) scalars, 50(3) 
scalars and fermions are 1, 2 and 3/2 respectively. The action takes the form 

S = -^T^J^ j dt tr {^D,X^DoXa - ^X,X, + ^[X„ X,]'..?j (2.1) 
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Due to the Myers term, the mass term and the commutator terms of 5*0 (3) scalars, the 
classical vacua are fuzzy spheres parametrized by Jj, the dimensional representation 
of SU{2) {[Ji,Jj] = i€ijkJk) , Xi = 2Ji, Xa = 0. There are a large number of vacua 
corresponding to partitions of integer N. In particular, there is a "trivial" vacuum 
Xi = 0, Xa = 0, which were conjectured from gauge theory analysis to be a single NS5 
brane wrapping 5"^ [9]. The gravity dual of this vacuum was studied in [10] (see also 
[18],[T9]), which is a linear dilaton background. The vacuum is dual to IIA little string 
theory on x R, providing further evidence for the identification of the single fivebrane. 
From the point of view of the M theory on asymptotically plane-wave spacetimes, these 
vacua are 1/2 BPS M2 and M5 branes preserving half of the supersymmetries of the 
asympotic plane- wave geometry. The gravity duals of these vacua were studied in [10] . 

There could be in principle tunneling between different fuzzy sphere vacua, exchang- 
ing some amount of DO branes, if we perturb the theory. From the 11 dimensional point 
of view, this process is the transferring of the longitudinal momenta of two fuzzy mem- 
branes. One naturally wonders whether there are instanton solutions in the plane wave 
matrix model interpolating different vacua. In fact, Yee and Yi [21] found a class of 
analytical instanton solutions interpolating between a fuzzy sphere vacuum Xj = 2 Jj and 
a trivial vacuum X^ = 0. There could be two approaches to the instanton action in the 
weak coupling gauge theory side. One is to study a bound in the Euclideanized action 
corresponding to the action of the instanton which was done in [21], and the other is to 
study the superpotential for each vacuum. The instanton equation can be read from the 
superpotential or by supersymmetry transformation condition^. 

We notice that in the weak coupling theory, the instanton solutions merely involve the 
three 5*0(3) scalars X^ and their superpartners, but not the fields in the 50(6) sector. 
The solutions for 50(6) scalars are Xa = for these instantons in the weak coupling the- 
ory. We can thereby write a "superpotential" [IS],[IZ] for this supersymmetric quantum 
mechanics, for each fuzzy sphere vacuum: 

W = ti (X.Xi + i^eijkXiXjXk) (2.2) 

It's easy to check that the lagrangian for X^ is correctly produced, and the action is 
bounded by the difference of the superpotentials between initial and final vacua: 



Sinst — — I dr ti 



9lj- 



oo 



9l 



/ + 00 1 
dr -[D^X.^dxWf -W\TJZ 



(2.3) 
(2.4) 



where we denote gymo/m? = g^, r = it. 

The instanton equation is the same as setting the square term to zero. 



D^Xi + dxW = DrXi + 2Xi + ie,,kXjXk = (2.5) 



^For example, the supersymmetry transformation conditions were utilized by f20' in classification of 
the BPS equations for the plane-wave matrix model. 
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The equation is the same as in [21] and a class of analytical solutions were found |21j : 



-^blocki^) ~ -^UpXripi''') ~ '^'^ripXrip ^ _|_ g2{r-To) C^'^) 



Here Xli^^^ is a rip x rip block matrix in X*, and we used = gauge. This solution 
corresponds to the transition where a fuzzy spheres with size at r = —oo gradually 
shrinks to zero and turns to a trivial vacuum Xli^^^, = at r = +oo. Moreover, we 
can have many such blocks labeled by p. So the instanton action is the difference of the 
superpotentials /\W between two vacua X*(±oo) = 2Jj|±oo, 



~ |(tr Jfl+oo-tr J2|_^) (2.7) 



9o ro ,=t^,3 3 

(2i 



1 np{nl - 1) 



where we used the identity that the second Casimir invariant Jf^^^ = Inxn- 

One can define the W for the trivial vacuum Xj = as some constant Wq as a 
reference point and then the superpotential for an arbitrary fuzzy sphere vacua can be 
defined as Wq plus the difference term in (12. 7p . Wq will be set to zero for convenience, 
because we only need to know the differences of W for different vacua. 

We can label each vacuum by A^2*'' copies of n'^'' dimensional irreducible representa- 
tions. These satisfy X]i^2*^^5*^ = N . In the large N^"^ limit, the instanton action will 
be simplified to 

In the supersymmetric quantum mechanics, the instanton amplitude between two 
vacua a and h is proportional to a factor given by the difference of two superpoten- 
tials e ^^^^t/n ^ ^ go _ \y"g set h = 1 in the second expression and also later 

discussions. This factor in the amplitude is largely suppressed if is small and the 
change Y2i |-^2*''^5*^^l-oo is large. However, there are fermionic zero modes around these 
instanton solutions, the path integral for the tunneling amplitude is zero, due to that 
the integration Yi J dCi of collective coordinates of fermionic zero modes C,i in the 

path integral gives zero. The vacuum energies for these vacua will not be corrected and 
they are exactly protected vacua quantum mechanically. Only fermionic operators with 
enough fermion fields which can cancel the integration of the fermion collective coor- 
dinates will receive instanton corrections to their vevs (see also [21]). In our case, the 
instantons preserve 8 supersymmetries, and there are at least 8 fermionic zero modes, 
corresponding to the massless Goldstone fermions of the broken supersymmetriesH 



^This is in contrast to the bosonic matrix models. If we had a bosonic model with the same action 
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The method for solving general instanton solutions in plane-wave matrix model and 
its moduli space has been discussed in detail by [23], [21]. The solutions for bosonic 
and fermionic zero modes around the instanton solutions have been analyzed by |21j . 
Both the instanton equation and their moduli space have similarities with that of the 
domain wall solutions in A/" = 1* SYM [26], [23]. The condition for allowed instanton 
solution between vacua X*(— oo) and X*(+oo) is determined by the initial and final 
representations of SU{2) and the gauge group. Specifically, it requires that the final 
representation modulo terms in the centralizer of it in the gauge group, falls into the 
orbit of the initial representation under the gauge group action, so that gX''{—oo)g~^ 
and X'^{+oo) are interpolated (fore more details see [23]). In the representation of Young 
tableaux, the condition for allowed instanton solution is that when comparing the initial 
Young tableau and the final Young tableau, there are some boxes moved from left to 
right, but there are no box allowed to move from right to left. An example is illustrated 
in figure [H The condition also implies that any fuzzy sphere vacuum can be interpolated 
with the trivial vacuum. 









N2 



N2- 1 



(a) (b) 

Figure 1: Illustration of the condition for allowed instanton solutions in the Young tableau 
representations: Comparing the initial Young tableau (a) and the final Young tableau 
(b), there are only boxes moved from left to right, but no boxes moved from right to left. 
In the example in the figure, a column with boxes moved to the additional boxes on 
the right of the first row. It corresponds to a instanton interpolating a fuzzy sphere with 
size A^5 to a trivial representation. 

Now we turn to the 2+1 dimensional super- Yang-Mills on 5*^ x i? , which is a limit 
from the plane wave matrix model. The 2+1 dimensional SYM on x R we discuss here 
is the theory when we reduce the A/" = 4 SYM on 5*^ x i? on the Hopf fiber of the |10j . 
The gauge fields in the A/" = 4 SYM reduce to the gauge fields in the 2+1 d SYM plus 
an additional scalar $. The S0{6) scalars Xa, a = 4, 5, 9 reduce to their components 
invariant on the S^. There are totally seven scalars. The gravity duals of its vacua were 
studied in [TU]. The consistent truncation of Kaluza-Klein modes from A/" = 4 SYM to 

of the SO(3) scalars but with no fermions, then the energies of each fuzzy sphere vacuum would get 
corrected by these instanton effects. The vacuum tunneling in a similar bosonic matrix model was 
analyzed in detail by 



6 



this theory was studied in detail by [27] 

The theory can also be obtained from expanding plane-wave matrix model around its 
fuzzy sphere vacuum, so under the continuum limit (e.g. [32]) of the matrix regularization 
of the large fuzzy sphere, it becomes a theory on S"^ x R. Suppose we start from the 
vacuum with A''2 fuzzy spheres, each of size ^ 1, we have N2N5 = N, where is 
the number of original DO branes. We expand around this vacuum. Fluctuations of Xi 
about 2Ji {N5 X N5) are decomposed into gauge fields on the fuzzy sphere and a scalar 
$ [9j|. The scalar $ describes the fluctuations of the sizes of the fuzzy spheres around N^. 

Due to its relation with the plane-wave matrix model, the vacua and their superpo- 
tentials analyzed above in the plane wave matrix model carry over to the U{N2) 2+1 
d SYM. For the U{2) theory, the vacua are quite simple, which are characterized by 
one integer n, this is the eigenvalue of $ ~diag(n, —n). The tunneling issue was briefly 
discussed in |9]. It was further studied in detail by [22], who presented the instanton 
equation and solved a class of analytical solutions for the U (2) theory. We will analyze 
the general U{N2) case, for we already know the superpontial for each vacuum via the 
relation to plane-wave matrix model. One can map the instanton solutions in plane-wave 
matrix model to those of the 2+1 d SYM theory. 

The action of the 2+1 d SYM can be written in two ways. One is in terms of scalar 
$ and gauge field A. Another is in terms of three scalar s Yi which are combinations of 
the $ and gauge fields A 

Y' = + e'^'^e^Ak, A, = e'^%e^ , $ = e'Y' (2.10) 

where Cj is a unit vector in i?^, ejCj = 1. The formalism in terms $ and A is more direct 
from a gauge theory point of view. The Yi originates in the plane-wave matrix model, 
from the fluctuations of the 5*0(3) scalars around the fuzzy sphere Yi = Xi — 2Jj. The 
formalism in terms of Yi is more convenient for seeing its origin from the plane-wave 
matrix model. 

We will write the action as in eqn (2.9) in [TU], further, we write it in such a way that 
it depends only on a dimensionless coupling constant g'^rn2/^i where m was originally the 
mass of 5*0(6) scalars and can be set to 1 for convenience. We rescale all fields so that 
the masses of 5*0(6) scalars, $ and fermions are 1, 2 and 3/2 respectively. The action 
takes the form 

S = j dt^d^n tr {^DoXaDoXa - ^X„X, + X,]^...^ (2.11) 

This is the same in effect as we set the radius of 5*^, l//i = 1/2 in [TO] . 

Because the vacua satisfy an equation / + 2$ = 0, where / is the 5*^ component 
of the gauge field strength and is quantized into integers time a half, the vacuum is 

•^The consistent truncation from M ~ A SYM to plane- wave matrix model were studied in for example 
[28] , [29] , [30] , [To] , where a phenomenon of the redefinition of the gauge coupling constant was observed, 
in both weak and strong coupling regimes. See also related earlier work in [5T] . 
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characterized by the eigenvalues of $ 



$ = idiag(4^\4^\...,nf^)) (2.12) 

after an unitary transformation [TO^. Notice that n"^^ = N^^^ — A^s {n^ <^ N^) are the 
fluctuations of the sizes of N2 fuzzy spheres. 

We want to use the superpotential in the plane-wave matrix model to define the one 
in this theory. Some of the mappings of fields between the two theories are 

Y, = X,-2J, tr^,,^.-.]V./d^Sl. [J„]^L, L, ^ -u„,eA (2.13) 

where Cj is a unit vector in R^, ejCj = 1. The superpotential can be derived from the one 
in plane-wave matrix model (12.21) as 

W = N5j d^n tTN,>cNAy^r + i^e.jkY^^jYk + leijkYiLjYk) + Wj, (2.14) 

where the second term is a large constant term 

W-h = ^''n.nJ^ (2.15) 

i=l,2,3 

which corresponds to the superpotential for the vacuum with exact N2 numbers of fuzzy 
spheres of size A^5, and the first term is the difference of the superpoential with respect 
to the vacuum $ = diag(0, 0, 0). The instanton equation is 

DrYi + dy.W = DrYi + 2Yi + lei.kYjYk + 2ieijkYiL,Yk = (2.16) 

The solutions to these equations in principle can be lifted from those solutions in the 
plane- wave matrix model for Xi, so the moduli space of these solutions are in principle 
known. 

We want to have also the expression in the variables of $ and gauge field A. As 
discussed in appendix of [10], that the 2-1-1 d SYM on R x can be lifted to a 3-1-1 d 
SYM on R"^ X S"^, which can also be obtained from reduction of a 4-1-1 d SYM on R^ x 
on the Hopf fiber of S^. Interestingly, [22] noticed that the instanton equation in the 
Euclideanized 2+1 d SYM on RxS^ is the same as the self-dual Yang-Mills equation on 
R^ X S"^ 122]. If we lift the Euclideanized 2-1-1 d SYM on an X3 direction, we can define 
the four dimensional gauge field strength on i?^ x 5*^ 

= fd'^n + {-2^)dTdx3, *J^={-2^)d^n + fdTdx3 (2.17) 

Note that the self-dual equation JF= * is equivalent to the equation of motion / -|- 2$ = 
0. The part of the action that does not involve the S0{6) scalars and is inherited from 
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the 50(3) sector in the plane-wave matrix model can be written as a | J JFA*JF term. So 
the action is bounded by 

Se = ^y"j^A*j^=^y"(j^-*j^) A(*j^-j^)+^y"j^Aj^ (2.18) 

^ \ j ^ ^ j d^^^^ (2-19) 

The bound is satisfied only when T = *T . Thus after matching parameters, we get that 
the instanton action should be 

Si^st = -\aW = -\l6N2N5{tT $'|+oo - tr $2|_oo) = - A^siVg V np|l^ (2.20) 
% 9o 9o 



This expression is of course consistent with the U{2) case studied by [9], 

From the 2+1 d SYM point of view, is absorbed into the definition of the coupling 
constant via the relation = El. So we have the instanton action for the 2+1 d 

9ym2 9ym0 

SYM 

S^nst = -^4iV2(tr $'1+00 -tr <l>2|_^) = -^iVaJ^npl^^ (2.21) 

92 ^92 j 

where (^1 = 9ym2/^y purely in terms of parameters in the 2+1 d SYM. 



3 Gravity description 
3.1 Euclidean branes 

In this section we turn to the gravity analysis of the instantons and superpotentials for 
the vacua of the two theories. In some regimes, the instanton solutions we studied in the 
previous section can be approximated as Euclidean branes wrapping some cycles in the 
dual backgrounds [TD]. 

The backgrounds dual to the vacua of these theories have an 50(3) and 50(6) sym- 
metry and thereby contain an and factors, and the remaining three coordinates 
are time, p and r], where p is a radial coordinate. These backgrounds are regular because 
the dual theories have mass gaps. The gravity equations of motion reduce to a three 
dimensional Laplace equation for V, it is in the space of p, rj if we combine p with an 
5^ [7]. The regularity condition requires that the location where the shrinks are disks 
(or lines if without the S^) at constant rj^ in the p, t] space, while the location where the 

shrinks are the segment of the p = line between nearby two disks [10] . V is regarded 
as a electric potential and there are charges on the disks. Due to fiux quantizations, the 
charge and rj distance are quantized. The charges are distributed in such a way that the 
disks are equipotential surfaces, and the charge densities on the edge of the disks vanish. 
The cycle we will embed the Euclidean D2-brane is the 3-cycle at p = between two 
disks, where the p = line combine the 5^ to form a 3-cycle. 
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(a) 



^ 



(b) 



Figure 2: The configurations before and after tlie transition. In (a) there is a large disk 
with N2 units of charge above the infinite disk. In (b) the large disk has A^2 ~ 1 units of 
charge and there is another small disk near origin with units of charge. There is a S3 
cycle that is along p = and between the disks, formed together with the S'^. Below the 
infinite disk, there are image disks. 

Since we use an Euclidean brane in a geometric background, such analysis is only 
an approximation when the transition is very slight or perturbative in nature. The 
geometries dual to two different vacua are different, so a non-perturbative analysis is 
generically needed, and it will be discussed in next section. The geometries before and 
after transition are approximated as the same background in this section. The wrapping 
of the brane should satisfy DBI-WZ action and global conditions. 

We first look at some regimes which obviously can be compared from both sides. 
We first look at the plane-wave matrix model. We look at in the gauge theory side 
transition from a vacuum with N2 copies of fuzzy spheres with size N^, to the vacuum 
with N2 — I copies of fuzzy spheres with size plus copies of size 1 (or a single 
NS5 brane with DO branes). We keep moderately large and N2 large. In the 
dual gravity description, each fuzzy sphere is a disk with charge proportional to A^25 and 
located at a distance r/g = f ^5 above the = plane (for more details, see section 2.2 
of [lO]). In the gravity picture, it corresponds to transition from a large disk with charge 
Q = \N2 above the rj = plane, to a large disk with charge Q' = \{N2 — 1) plus a very 
small disk near the origin (see figure [2]). We consider the small disk and its image as a 
small dipole with dipole moment 2{Q — Q')^q = near the origin, and that they do 
not back- react to the geometry. The total dipole moment, which corresponds to N DO 
branes, is conserved. 

There is a nontrivial 3-cycle at p = 0, between the large disk and the rj = plane. 
The cycle has topology of an 5*^, with the S'^ shrinks only at r/ = and f] = fJo- The 
Euclidean D2 brane wraps the S"^ and its Euclidean time r is embedded along rj direction 
extending from r) = to rj = rj^. We should find r]{T) as a function of r from the DBI- 
WZ action. Such Euclidean branes preserve half of the total 16 supercharges, since they 
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extend along rj and have an additional projection condition for the Killing spinors. This 
Euclidean D2-brane has Euclidean DO-brane charges, and it mediates the transferring 
of one unit of D2 charge and units of DO charges, see figure [2l We utilize the general 
background solution in eqn (2.20-2.24) of [10]. It is parametrized by an electric potential 
y as a function of p, rj. 

The DBI part of the action is 

SDBi/r2 = - j dTd^Qe-"^ ^det ^/det G± + det(27ra'F - B) (3.1) 



2(3.2) 



where T2 is the Euclidean D2 brane tension, g^, g^jr^, g22 are various metric components on 
the t, 77, and 5*^ directions in the background fTO]. Gy, Gx_ are the time component and 
the spherical components of the pull-back metric, and det = 922 2TTa'F — B = 

{irN^a' — 2rj— ^^' ) smOdOdcj) according to flux quantization. We have chosen the gauge 
that 2TTa'F — B is zero where rj is at the top disk. Dot and prime are the derivatives 
w.r.t. logp and rj. The WZ part of the action is 

Swz/t2 = - J[C3 + {27Ta'F-B)ACi] (3.3) 

r AV • • ' 1 

= - drAn-. -[-V - V (-TiN^a' - r])] (3.4) 

J V-2V 2 

Due to the background flux, the Euclidean D2 brane worldvolume gauge field has 
a source of units of magnetic charges, and these are precisely cancelled by the DO 
branes ending on it. This is consistent for wrapping the brane on this cycle. This is 
similar to the situation of Euclidean D2 branes wrapping on the SU (2) group manifold 
with flux discussed in [33] . 

The expression of the action can be simplified if we expand the electric potential near 
p ~ 0, as 

^= E (-l)^"|J^>- = -A-o + fp'-|^P' + -. A-.A-J (3.5) 

i=0,even ^ ' 

where K is only a function of rj. K has some properties such as > 0, > 0, in 
order to satisfy the regularity condition for the geometry. The action is then 



47rr2 



-jdr ( ) ^ \J\K + K'i^N, - r^)Y - \kK"{7.n, -Vy]l^ + {d^vY 
dr=^[K + K'{^N,-v)] (3.6) 
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We have set a' = 1. Since the action does not exphcitly depend on r, we can have a 
conserved quantity H = PdrTj — L = const. = E, where P = dQ^^jL. It is obvious that 
the constant solution rj = is a. special case of a class of solutions to the equations, 
and from this solution we get E = Q. Setting = 0, we get the equation of motion 

2m_N, - ri) 

The solution for the Euclidean D2 brane is 

= (3.8) 

There could be an overall constant on the RHS, but it can be absorbed by shifting r. 

We plug in the solution back into the action and integrate out to get the final answer 
of the Euclidean D2-brane action 

Se = -4^^2^' ' 2K{'^N,-r])dr] (3.9) 

= --[V{Vo)-V{0)-VoV'm (3.10) 

TT 

where we used integration by parts several times to arrive at the last expression, and we 
used T2 = a' = 1, tiq = \N^- Vijj) in the expression is understood as the potential 
V{p, 1]) along the p = line. The expression (13.101) has the right property of the invariance 
under the shift of \^ by a linear term in rj, since the gravity solution should be invariant 
under this shift. 

To get more intuition, let's look at the example of two infinite disks, which is the NS5 
brane soultion in ^U\. We wrap a Euclidean D2 brane on the throat region of the NS5 
brane solution. For that solution |10j . 

V^ = i/o(r)sin^, r = |f, 6 = ^ (3.11) 

where go is a constant, so = —V"\p=o = -^gQ^sinO. The action of the Euclidean D2 
brane is 

Se 



J dr^/{n - 0)2(c2 + 1) + 4sc(7r -9) + 2s^^/8 + (9,0)2+ J rfr4[s+c(7r-0)] 

(3.12) 



SirgQ V2 

where c = cos 6', s = sin 6, 9 G [0, vr], and the equation of motion is 



(3.13) 

s + c{tt -d) ^ ' 



The solution for the Euclidean D2-brane is then 

-2r _ -2roSinfvr 



e"^^ = e-^^« — ^ — (3.14) 
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where tq is a constant. When 6 = it, r(7r) = tq, when ^ = 0, r(0) = +00, so this tran- 
sition corresponds to transferring a small charge (corresponding to one unit of A^2) from 
the top disk to the bottom disk. The Euclidean D2 brane action for this solution is 

= To- 

Of course, the solution in the form (13.101) works for general configurations of disks 
in the plane-wave matrix model, when there are nearby two disks. The Euclidean D2 
branes wrapping on that cycle mediate transferring of charges between two disks while 
keeping the total dipole moment fixed. It also works for the 2+1 d SYM. In that case, 
the Euclidean D2 brane mediates similar transferring of charges between two nearby 
disks, while the total dipole moment of the system is kept to zero, because the theory is 
expanded from the plane-wave matrix model and the disks are already in the center of 
charge frame. 

There are issues of the range of validity that need to be addressed. One is that 
the calculation of the embedding of the Euclidean D-brane needs that the gravity back- 
grounds are weakly curved, this requires that the two disks where the cycle is in between 
are quite large, and their separation is also quite large. This needs the parameters 
i.9ymo/^^)^2, ^5 to be relatively large. On the other hand, had better not be too 
large, this is because from the Euclidean D2 brane point of view, there are N^, DO branes 
ending on its worldvolume, if Nr, is too large we would consider these DO branes 
as a single NS5 brane, and the brane configuration should be changed. From the disk 
point of view the small disk near the origin should be much smaller then the top disk, 
so A^5 ^ And these conditions have already guaranteed that the charges on the 
Euclidean D2 brane are much smaller than the background charges. 

A second issue of the range of validity is that the Euclidean brane calculations only 
describe the transition between very close vacua or very close backgrounds. The geometric 
backgrounds before and after the transition are approximated as the same geometry. So 
in principle we need non-perturbative description which could interpolate between very 
different backgrounds, as long as it is allowed by the condition in section [2j This will be 
discussed more in section 13.21 in terms of the non-perturbative superpotentials. 

The condition for the instanton solution as discussed in section [2] in the gauge theory 
side is that there are boxes moving from left to right only, in the Young tableaux. The 
number of boxes are the number of DO branes transferred. This condition in the gravity 
side in our case is the transferring of charges from the higher disk to the lower disk. The 
process in figure [2] corresponds to the illustration in figure [H 

The instantons and Euclidean branes we studied so far describe tunnelings between 
spherical giant gravitons and point giant gravitons in the lid plane- wave geometry. This 
can be related to the limit of those instantons tunneling between spherical giant gravitons 
and point giant gravitons in AdS^^j x S'^''^ studied by [35], [36], (see also [37]). The 
difference is that the giant gravitons ([34j,[35j,[35]) in their analysis of the tunneling 
are probe branes in the AdS4j x S"^'^ geometries and in our case we have replaced the 
giant gravitons with the back-reacted geometries taking into account the back-reaction 
of them on the plane-wave geometry. It is analyzed both in [53] and [2T] that these 
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instantons have 8 supersymmetries, and are 1/4 BPS objects from the point of view of 
the asymptotic geometries. In our case, the kappa symmetry condition on the Euchdean 
D2 brane will impose an additional projection condition, with respect to the one for the 
backgrouncEl- Again, there will be at least 8 fermionic zero modes corresponding to the 
broken supersymmetries and the fermionic zero modes suppress the tunneling amplitude 
and forbid the mixing of the vacua. 

One can also consider an Euclidean NS5 brane wrapping a six-cycle enclosing the 
disk, which is the Eg cycle in [10] . This Euclidean NS5 brane has N2 units of Euclidean 
DO charges. It interpolates between a vacuum of iVs NS5 branes with size N2 and the 
vacuum with A^s — 1 NS5 branes with size A^2, plus A^2 DO branes. On the gravity side, 
the transition corresponds to lowering the disk by one unit of the height. This instanton 
solution is allowed from gauge theory side but its action is zero in the weak coupling limit. 
From the gravity side and strong coupling point of view, the action should be non-zero. 
The existence of this Euclidean NS5 brane is a quantum effect from gauge theory point 
of view. Since the upper surface and lower surface of the disk are treated differently, it 
is more convenient to use the coordinates of the Xi = Xu, X2 = V, y = V appeared in 
the Toda equation [7], [10]. From lid point of view, it's a Euclidean M5 brane wrapping 
5*^, with its worldvolume r extending along X2 of the M5 strip at y = and traveling 
along the xi circle, with N2 units of momenta Px^ conjugate to xi. From the point of 
view of the IIA Euclidean NS5 brane, there is a worldvolume scalar corresponding to xi, 
and the action [38] can be written aaj 



T = {drXi-Ct)dT (3.15) 

where Ts = (2^5 ct' = 9tt, 9x2X21 955 ^^^^ the components of the metric along the t, 
X2, directions, -86,^*5 are the pull-backs of the potentials dual to -B2, C*3, and JF is 
the worldvolume one-form flux. We look at solutions that P^, = ^^^^ ^ = that is 

■^1 d(aTXi) Rii ' 

quantized into integers A'2 corresponding to DO charges, and then the equation of motion 
will yield X2{t) as a function of r, which characterizes the embedding of the brane. If 
we expand the solution near an isolated disk, we can approximate it by the solution of a 
single disk, which was the solution (2.49)-(2.54) of [10] corresponding to the 2+1 d SYM. 
Using this solution, we then estimate the action to be 5"^ ~ 9ym2{9ym2^'iY^^ ~ 9yTnoVoPoi 
where //q and Po are the height and size of the disk, and gymo = go is the coupling 



"'in this case, we have the kappa symmetry projector on the Euchdean D2 brane: Te = 
{^jm^^iTrTgT^ - ^rgr^rii)e = e. Here ^^22,322 are the S'^ components of the worldvolume 2- 

form flux and the pull-back metric, = EjT,^ = Lr^drTlTjj+iVj) / .flT^Hd^ , e = S{x'-)eo, 



-gtt ' ' 'I ^" \l -gtt 
where eq is a constant spinor subject to projection conditions r-jFueo = iTjjT-gT-^eQ — eo- 

^Note if the NS5 brane with N2 DO charges wraps instead of , the configuration can be T- 
dualized to {1,N2) 5-brane in type IIB theory after taking T dualities five times [3^. The IIA NS5 
brane action with D2 or D4 charges have been studied in detail in for example [40] , |41J . 
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constant in the plane-wave matrix model. This is a quantum effect and strong coupling 
result which is not seen in the weak-coupling gauge theory. 

3.2 Electric charge system and superpotential 

Now we analyze what the quantity W is in the electrostatic picture. In the gauge theory, 
the instanton action is proportional to the difference of the superpotentials of the two 
vacua as defined by Smst = —^'^W in (12. 7p . In the case of the two vacua that is very 
similar to each other, we used the approximation of an Euclidean brane action to match 
the instanton action. Of course, the Euclidean brane action will be the strong-coupling 
result. For the configuration of two vacua studied in section [STTl involving the exchange 
of a unit of charge (corresponding to a unit of N2) between two disks, we arrived at the 
expression 

Se = --[Vivo) - ViO) - VoV'm (3.16) 

71 

This quantity is proportional to the change of the total energy U in the electrostatic 
system (including image charges). Note that the relation between D2 brane number and 
electric charge is Q = [lOj, so the charges are in units of the charge transfer in 

2 

this process is AQ = The energy that the unit charge and its image charge released 
is 

- Af/eharge = 2AQ[V{rio) - V (0)] (3.17) 

On the other hand, the total dipole of the charge system is conserved, the dipole had 
been lost by the amount 2AQriQ by transferring the unit charge from the top disk. After 
the transferring process, we have the small disk (and also its image disk) near the origin, 
as in figure [2](b). We make an approximation when a pair of small disks near the origin 
is approximated as a point. We thereby have a small dipole with dipole moment 2AQriQ 
at the origin. The electric field at the origin is —^'(0), which is pointing upward, so the 
coupling between the dipole and the electric field there gives an negative energy, this is 
the energy the dipole released: 

- Af/dipoic = {2AQ7]^)[-V'm (3.18) 

So the total energy difference of the system between final and initial configuration is 

Af/totai = -2Ag[V(r/o) - V(0) - r/oV'(0)] (3.19) 

which is just proportional to the Euclidean D2 brane action, in appropriate unit. Of 
course this definition of the change of the total energy is invariant under the shift of the 
potential by a linear term hrj. The energy difference will not be changed due to this 
shift of V . This is because the charge AQ, will release an additional energy 2AQhriQ, on 
the other hand, due to that there is an additional electric field downward with magnitude 
6, the dipole's energy is increased by 2 AQr/gfe, which exactly cancels the additional energy 
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released from the charge. So this expression is consistent with the definition of the change 
of the total energy. 

We therefore can identify (in the strong coupling regime) Smst = Se = ~3^U, and 
thereby the superpotential for each vacuum can be defined in terms of the energy of the 
electric configuration corresponding to that vacuum 

U = -^,W (3.20) 

up to an overall constant shift which is not important when comparing differences. The 
minus sign is due to our conventions. So the superpotential of the system is proportional 
to the total energy of charges (including image charges) 




(3.21) 



i 



and Wo which corresponds to the single NS5 brane vacuum is set to zero for convenience. 
i represents different groups of charges at t] = rj^^ with potential Vi. Now we see that 
the factor of the instanton amplitude is given by the energy difference between the two 
configurations of the electric charge system or the eigenvalue system. The factor in 
the instanton amplitude is largely suppressed if the energy difference is large. Note 
that there are still other parameters like the a' and h, which have been set to 1 for 
convenience. These charges are associated with eigenvalues from the S0{6) scalars and 
will be discussed in section HI 

It might be useful to notice that the expressions 03.211) can be reduced to integrals of 
forms in p, r] space. We define one forms locally dV = dr^Vdr] + dpVdp, dq = 2np *2 dV = 
2'K{—dr^Vpdp + dpVpdf]), where *2 is the fiat space epsilon symbol in p,!] space. The 
Laplace equation for V is just d{dq) = 0. 

We choose one-cycle Ilj enclosing each disk, and one-cycle Fj between each disk and 
a reference disk (whose potential is set as zero). We see that 

Q^= [ dq, V,= f dV (3.22) 

Since Ilj cycle and the forms a 6-cycle, and Fj cycle and forms a 3-cycle in 9 spatial 
dimensions, the Ilj and Fj cycles do not intersect, Ilj fl Fj = 0. We can thereby form a 
two-cycle by the cup product of the two cycles, so (13.211) can be written as 

W = -^[ dqAdV = -^ [ \VV\''2TTpdpAd7] (3.23) 

where W is the gradient of V in p,ri space. So we have defined the superpotential of a 
boundary gauge theory by purely forms appeared in the bulk of gravity duals. This view 
is very similar in spirit to the definition of the non-perturbative superpotential in 3-1-1 d 
J\f = 1 gauge theories via a fiux background after geometric transition, e.g. ^2] . 
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Let's compare the superpotential in the weak couphng case. The charges in the 
electric system originate from the eigenvalues of the matrices in the 5*0(6) sector in the 
gauge theories. We will discuss this more in section HI Classically, these eigenvalues do 
not interact with each other. They distribute under the external electric potential, which 
are quadratic (~ p^) in the gravity dual of either theories [lO]. So they are sitting on top 
of each other at fixed position 77^ at p = 0. Quantum mechanically, the eigenvalues repel 
each other and they are not coincident. They thus form extended disks. The classical 
and quantum pictures, correspond to weak coupling and strong coupling effects for the 
superpotential. In the weak coupling analysis, the superpotential does not involve any 
dependence on the 5*0(6) scalars, so this must get corrections when coupling increases 
and becomes larg^. 

For the plane- wave matrix model, the background potential is Vf, = —{p^Tj — \rf)^ 
where Qs = in'^gyj^o (^-g- [IS]) and gs is the string coupling that appears at the asymptotic 
DO brane near horizon geometry |31]. Since in the gravity side, m is the unit of energy 
for [To], which is 1 in our case, so we have m = 1. In the weak-coupling regime of the 
gauge theory, since the charges sit at p = and r] = rj^, they feel potential Vi = — ^^igi vh 
so the superpotential from (13.211) is: 

w=^.j:Q^v!=j:l^^''Nt (3.24) 

i i 

where we used = N^^^ = ^ [10] from flux quantization conditions. This is the 
same as the field theory calculation in weak-coupling regime, cf. (12.91) . The factor of 
the instanton amplitude is largely suppressed if gs is small and the change of the energy 
of the system is large. 

The vacua of the 2+1 d SYM are inherited from that of plane-wave matrix model. In 
the disk configuration of the plane-wave matrix model, if we look at the geometry locally 
near a group of disks that is isolated from all other disks, including their images, then the 
geometry locally approaches to the solution dual to 2+1 d SYM theory. The vacuum of 
the 2+1 d SYM thereby is also characterized by the total energy of the electrostatic sys- 
tem. A special case is that when we start from plane-wave matrix model and look at the 
configuration where all the disks are very high above 77 = plane and they are relatively 
near each other. Their average position is t]^ = f A^s. We define new coordinate relative 
to the average position rj = 1] — tiq = \n^, (jj <^ ?7q) and = N^ — N5, [n^ <^ N^). The 

total number of the D2 brane charges is = X]i=disfcs -^2*^- coordinate 
frame, this is the configuration dual to the U{N2) 2+1 d SYM. is proportional to the 
location of each disk in the new frame. The case with all ns = corresponds to a single 
disk at r/ = 0, and is the simplest vacuum with unbroken ?7(A^2)-The number is used 

^Some proposals for the quantum mechanically corrected action of the plane- wave matrix model 
taking into account the quantum effects of the SO{'o) fields have been put forward, for example, in 
[45], [46] using fuzzy five-sphere constructions. It is also possible to use the approach of [47] . 
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(a) 



(b) 



Figure 3: Comparison of strong coupling and weak coupling superpotentials. In (a) the 
charges or eigenvalues are extended into disks due to eigenvalue interactions in strong 
coupling regime. In (b) they are on top of each other at p = line and not extended 
in the S0{6) directions in the weak coupling limit. The superpotential is given by the 
energy of the electric charge system. 



in the matrix regularization for the S*^. The electric potential near this group of disks is 
V = 4^{p'v - Iv') = - 2v') - + - (3.25) 



They basically feel the same external potential ^^p-(p^ — 2?]^) as in the gravity dual of the 
2+1 d SYM theory discussed in [10], with a large constant shift due to our expansion. The 
dot terms are the linear term in rj and higher order terms in rj, which are not important. 

In the weak coupling regime, the charges are not extended in the p direction. The 
total energy of the system (excluding the large constant term) gives the superpotential 
from the relation fl3.2ip 

= E «■ = E (3.26) 

The large constant term from the 7]^ term in the electric potential, corresponds to Wj- = 
I tT^NxNJf term for the fuzzy sphere we are expanding around in the plane- wave 
matrix model, cf. fl2.15p . 

We thereby find that the weak-coupling superpotential analyzed in section [2], which 
does not involve 5*0(6) scalars, is the situation in the gravity side when the eigenvalues are 
not extended in the p direction, or the 5*0(6) directions. This is when the eigenvalues do 
not have interactions between each other. To bring in the interactions, we need quantum 
effect, and this is shown in the strong coupling regime in the disk picture. These are 
illustrated in figure [3l 

Now we briefly discuss some technical aspects of solving the value of the superpotential 
in the strong-coupling regime, for the disk configurations. We first discuss the situation 
with a single disk above 1] = plane (13.101) . and then the situation with some arbitrary 
numbers of disks (13.211) . 
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We first discuss the situation of a pair of positively and negatively charged disks 
under the background potential Vb = ^(p^?7 — |?7^)- The size of the disks is Po 
their separation is 2r)Q. Apart from the background, the charges have mutual repulsions 
from other charges on the disk and attractions from the image disk. Suppose the charge 
density and the potential on the disk is a{p) and Vq. We can scale out everything and 
the problem then is only determined by the ratio k, = y-. We make the rescaling y = 
V/Po, X = p/po, a{p) = j-^plg{x), V = j-^plv. 

The integral equation in these variables is (which involves dimensionless functions 
only) 

,,9 2 o, f'^'^ f^, q(x')x'dx'd(f) q{x')x'dx'd(b 

{{x^K — K^)+ / / [ ^ ' ^ ' =]}v^g[o,i] = const. 

3 Jo io V — 2xx' cos + a/x^ — 2xx' cos + x'^ + Ak^ 

(3.27) 

with the supplementary condition that (7(1) = due to vanishing of the charge at the 
edge of the disk. When g{x) is solved, it is also a function of k. 

The integral equation is very difficult to solve due to its complicated kernel. Provided 
we know the function g{x), since we already know the total charge of the disk, what 
we need to compute is the potential V on the disk. Since the disk is an equipotential 
surface, we can know it from the potential at p = 0, this consists of a contribution from 
the background and a contribution from the charges. The value of the superpotential 
dMH) is 

W/{^-N^Nl) = 1 - ^rtTT ^ dx g{x 



(3.28) 



(s2 + 4/t2)i/2_ 

It contains two terms. The first term "1" on the RHS gives a contribution to the su- 
perpotential that is the same as the weak coupling gauge theory result W = ^N2N^, 
it's the result if the charges were all sitting at p = 0. The second contribution is from 
the integral term, resulted from the effect of mutual interactions among charges. We do 
not assume that k has to be small, so in principle, if we were looking for perturbative 
corrections to the weak coupling answer, the integral term gives the corrections relative 
to the first piece, in the large k regime. The k depends on N2 in a very non-linear 
way. In the small k limit, it is estimated in [10] that pg = ^n'^gy^QN2, rj^ = jN^, so 

The problem of solving the integral equation has been simplified by Abel transforming 
the kernel in (13.271) as done in pIB]- The new integral equation is 

fix) - f - ,2^( ,,2 f(t)dt = 1 - 2«^' (3-29) 
J_i n Ak^ + [x — t)^ 

where f{t) is the Abel transform of g{x), f{t) = ^ ^T^^^)^' ~ "^o + f^^^, with the 

additional condition /(I) = 0, f'{l) is bounded [18], and a is a constant determined 
by these boundary conditions. The new kernel is the one appears in the Love's integral 
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equation [H] , [19] , [50] . The solution in this form seems very promising. The kernel can 
even be simplified a little bit more in the small k limit [49j . 

Now let's discuss the calculation of the superpotential for the electrostatic configu- 
ration with an arbitrary number of disks. Once we specify the charges Qi of the disks 
and their heights 77 j, what we need to know is the electric potential of each disk Vi, in 
order to compute the superpotential (13.211) . We can thereby use the conformal transform 
techniques in [TO] . 

In [To], the complex variables were introduced: w = 2dzV = dpV — idrjV and z = 
p — Po + if]. In the large disk limit, the general solution for the plane- wave matrix model 
is 

^-^ = fl l'^~'%^ H»»+i) (3-30) 

where there are n finite disks located at iaj, and n fivebrane throats located at zcj, and 
finally the infinite disk at ian^i. We can integrate it to get the potential of each disk 

wdz = —Re / IT 7 ^(—ian+i)wdw (3.31) 

and then plug this in (13.211) . The rest is to figure out the relation between ai,Cj and 
iV«,iV« [TO]. 



4 Emergence of gravity picture and eigenvalue sys- 
tem 

In the previous sections, we find that the instanton action is associated with the super- 
potential of each vacuum, which is in turn related to the emergence of the geometries. 
The geometries are emergent from matrices in these theories. There are different types 
of emergence in our cases. The emergence of odd dimensional spheres and the even di- 
mensional spheres seem to be different in some cases. In this section, we will discuss 
qualitatively the emergence of the extra coordinates and the electrostatic picture. 

For the plane-wave matrix model, we have to explain the emergence of 9 directions, 
5*^, 5^, ?7 and p. We first discuss the emergence of S*^ in the plane- wave matrix model. 
This is by the SU{2) commutation relations of Ji (e.g. t32j). Other higher even dimen- 
sional spheres may similarly emerge in theories constructed using higher even dimensional 
fuzzy spheres for example [5T].[52j. 

These fuzzy two-spheres have radial motions, which are characterized by the scalar 
$, most conveniently seen in the 2-|-l d SYM. Comparing sections [2] and [31 we find that 
the eigenvalues of $ are proportional to n^^\nf \ ...,nl^^\ which map to the heights of 
each disk in rj direction: ri^,ri2, f]^^. So the 77 direction is emerged from the eigenvalues 
of $. 
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Next we come to the emergence of the odd dimensional sphere and the direction 
p. From intuition, due to the reduction from A/" = 4 SYM, the SO {6) sector in these 
theories are quite similar to A/" = 4 SYM. So the emergence of should be similar to 
that of A/" = 4 SYM 

In order to address this problem, we embed these vacuum geometries into the larger 
sector of 1/8 BPS geometries. These 1/8 BPS geometries are dual to 1/8 BPS states in 
these gauge theories. We can look at 1/8 BPS states satisfying E — {J a + Jb + Jc) = 0, 
where E is the energy and Ja, Jb, Jc are the three U{1) i?-charges associated with the 
excitations by three complex scalars from the 5*0(6) scalars. These 1/8 BPS states are 
in turn embedded in an even larger Hilbert space. Finally our BPS vacuum geometries 
are the ground states in the 1/8 BPS sector. 

Though quite similar, a difference of the theories we analyzed, relative to AT = 4 
SYM on X R, is that we have many ground state geometries characterized by a broken 
gauge symmetry U{N) Y[U{N^'^'>). The most convenient theory is the 2+1 d SYM 

i 

where a particular vacuum is specified by the eigenvalues of the scalar $. If there are 
iV*^*^ coincident eigenvalues, then there is a [/(iV*^*)) symmetry. These eigenvalues form an 
extended disk in the electric charge picture, but this could have alternative descriptions 
if we embed the geometries into larger sectors. If the disks are far away from each other, 
the repulsion of charges on each disk is quite similar to the case of a single disk. In this 
limit, we study only the case of a single disk, that is, we look at the block of the SO (6) 
scalars corresponding to A^*^*^ and only excite U{N^'^'>) invariant 1/8 BPS states. Then 
finally we can excite all the blocks and build Y[ U{N^^'^) invariant 1/8 BPS states. 

i 

For our purposes, we can consider first a single disk vacuum in the 2+1 d SYM with 
unbroken U{N) symmetry, and then other cases will be quite similar. These 1/8 BPS 
states can be described by the harmonic oscillators of the 5*0(6) scalars and only s-wave 
modes of the 5*^ are relevant for this class of BPS excitations. We use the methods similar 
to those of [1] . [5] . [Sj^. The action that is relevant for these excitations of 5*0(6) scalars 
is the matrix quantum mechanics 

S = [ dt tr (iXaXa - IXaXa + ^ [X„ xA (4.1) 



m 



9ym2/ 

where a, 6 = 4, 5, 9 and we have set /i = 2 in the lagrangian. We set = and have 
a Gauss law constraint on U{N) gauge singlet physical states, ^ = 0. Furthermore, 
the commutator term^ will not be included in this class of BPS excitations because they 



^One can also use the approach of collective fields for the three complex matrices, for example in [53], 
although in the loop variables the other two matrices are treated as impurities with respect to the first 
matrix, the method is in principle general. 

^Strings can also emerge on the backgrounds of these emergent geometries, when we turn on the 
commutator terms. They are fluctuations of the impurity matrices on the background of other matrices 
(corresponding to the background geometry), see the approaches of, for example, [54] , [55] , [56] , [57] from 
gauge theory side. 
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contribute positively to the energy but not the i?-charges. Then we define the conjugate 
momenta Pa for each X^. The energy and -R-charges are 

E = hliPaPa + XaXa), Ja = tl {P,X, - P,X,) , (4.2) 

Jb = ti (PeXj-PrX,), Jc = tT{PsX,-P,Xs) (4.3) 

We can form the usual three complex scalars Z^, and their conjugate momenta Ul = 
—i^Y' (-^ ~ They satisfy the standard commutation relations [2'™", III'" " ] = 

ihSnm'Sn'm, = 0, SO wc cau define creation/annihilation operators 

ai = ^{Zl-^U[), b[ = ^{Zl + ^U[) (4.4) 

and their conjugates. The energy and i?-charge operators are then written as 

E = tr^(a[aL + 6[6L), Jl = ti {a{aL - blbi) (4.5) 

L 

So it's clear in order to get 1/8 BPS states, we should keep only the three oscillators but 
not the three b^j^ oscillators. We have then in this 1/8 BPS subspace of the Hilbert space 
a three dimensional harmonic oscillator Hilbert space. Its phase space is six dimensional. 
The singlet condition tells us that we should look at states of products of traces 

k 

\^)=U^r{a\y^{a^^r^{ayr\0) (4.6) 
i=i 

with lj,mj,nj ^ A^. Because the commutator terms vanish, by unitary transformations 
we may simultaneously put Za, Zb, Zq into the form of triangular matrices, each of which 
is a sum of a diagonal matrix and an off-diagonal triangular matrix. We integrate out the 
unitary matrices and also the off-diagonal triangular matrice^. Let the diagonal matrix 
elements be ZAi-, zbi, ^a- The wave function corresponding to (14. 6 p is of the form 

k 

j=\ i i i 

up to some normalization factors, including the factor from integrating out off-diagonal 
matrices. l\{zAi-,ZBi-,zci) is a Van-de-monde determinant from integrating out unitary 
matrices. See similar wave functions in [1], and its 1/2 BPS case [5] . [6] . [58] . [59] . [60] . 
There is an interesting property. The wave function, besides the universal Gaussian 
factor, is a holomorphic function in ZAi^ZBi^zci^ due to the 1/8 BPS condition. This 

^We refer the readers to [59], [30] for more detailed discussions and more subtleties for the case of one 
holomorphic complex matrix. 



22 



holomorphicity provides some simplifications in the wave functions. Wave functions 
that depends also on ZA_i,ZBi,zci are non-BPS. The complex structure is related to the 
symplectic structure in the phase space C^. Each zu coordinates in the phase space are 
paired by a physical real coordinate and its conjugate momentum, zn = xu + ipu 

In the large limit, the wave function has a thermodynamic interpretation, the 
wave function squared describes the probability of a particular eigenvalue distribution 
in phase space [Ij. The saddle point approximation of the wave function gives a droplet 
configuration in the phase space. The polynomial part of the wave function provides 
repulsion forces between eigenvalues. The Gaussian factor gives a spherically symmetric 
quadratic potential. The ground state should be a spherically symmetric droplet with 
5*0(6) symmetry in C^, which is bounded by an S^. Because of the spherical symmetry 
of the ground state droplet, the emerges in the gravity description. We see that the 
eigenvalues are extended in the radial direction of R^, the p coordinate in the gravity 
description is really mapped from the radial direction of the phase space R^. In the 
electric charge picture, charges are extended along p, if we tensor the p with the S^, 
the charges are really distributing in a six dimensional space. The electric disks emerge 
from the eigenvalues of the SO (6) matrices. This would be more clear if we had all the 
regular 1/8 BPS geometries of these theories. The then will be a particular enhanced 
symmetry for only the ground state geometries. 

If we start from 1/4 BPS sector or 1/2 BPS sector, we will have similar but reduced 
droplet descriptions. They are described by two or one dimensional harmonic oscillators 
similar to the discussion above. They will have a or phase space. The ground 
state droplet will have S0{3) symmetry in R^ bounded by an or U{1) symmetry in 
R^ bounded by an S^. In the electrostatic picture, there is an isometry from the S^, 
that could combine with p into a real two dimensional disk. If we had all the 1/2 BPS 
geometries for these theories, they should have a S*^ isometry out of 5*^. The ground 
states have an enhanced symmetry of which combines p to form a disk. Small ripples 
on the edge of the disk then describe BPS particles traveling along the equator of the 
5^ We have a further evidence from the 3+1 d U{N) A/" = 4 SYM on i? x S^/Z^ with 
k/N small. In the electrostatic picture, the vacuum are periodic disks. If we look at the 
vacuum with one disk in a single period, because k/N is small, they look like charges 
in a cylinder, and the Laplace equation will be two dimensional in the space of p and 
an S^. On the disk, the charges feel logarithmic repulsions plus an spherical symmetric 
quadratic potential. The integral equation for the equilibrium configuration is exactly 
the same as the one in eqn (2.16) of [4J for the 1/2 BPS ground state droplet of A/" = 4 
SYM on R X . We see that their descriptions are consistent and should be exactly the 
same for k = 1. 

These ideas should be compared with the gravity results. A class of 1/4 BPS and 
1/8 BPS geometries in jV = 4 SYM, relevant to the matrix model discussed above, 
have been studied by [HI], [62] and by [63]. They have U{1) x SO (A) x S0{2)r [61j,[62] 
and U{1) X S0{4) symmetries [53] (after analytical continuation of AdS^) respectively. 
A Kahler structure were observed in both cases. In the 1/4 BPS case, it was reduced 
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to a five dimensional Monge-Ampere type equation which is highly non-linear. The 
boundary condition is roughly to divide regions where the 5*^ or 5*^ shrinks smoothly. 
Thereby there could be a droplet configuration in a four dimensional space. However, 
the regularity condition and the topology of these solutions are very complicated. To 
understand the differential equation in terms of the emergent forces between eigenvalues 
is a very challenging problem. Our case of the emergence of the electric charge picture 
of eigenvalues from gauge theories may bear some similarities to these more complicated 
cases. 

The emergence of the odd dimensional sphere from phase space discussed above can 
be generalized to other cases, for example S'''. We consider the emergence of AdS4 x 5*^ 
from M2 brane theory on RxS^. Apart from the time and S"^ that are already present, we 
need to explain the radial coordinate r in 74^5*4 and the S'''. We can look at the 1/16 BPS 
states of the M2 brane theory, by looking at a BPS bound in which energy is bounded 
by the sum of 4 U{1) i?-charges out of 5*0(8)/^. They correspond to excitations of four 
oscillators associated with the 4 complex scalars for the 1/16 BPS chiral primary states. 
The 1/16 BPS sector has an eight dimensional phase space. The ground state droplet 
should have an enhanced symmetry of S*^ in the phase space. Similar picture for the 1/2 
BPS sector of M2 brane theory has been studied in [7], the emergence of the 5*^ (part of 
5^) for the ground state in a two-plane has been observed. 

Other ways of emergence of odd dimensional spheres, similar to the even dimensional 
ones, have been studied by for example [6l],[52],[65]. Their emergence is similar to the 
5*^ in our cases. 

5 Conclusion and discussions 

In this paper we studied the vacua of gauge theories and the dual gravity descriptions of 
each vacuum. These theories are U{N) plane-wave matrix model and 2+1 dimensional 
SYM on R X S"^. They have a large number of supersymmetric vacua. Each vacuum 
can be consistently lifted to the Af = 4 SYM on i? x 5*^, where they become pure gauge 
transformations of the vacuum of the A/" = 4 SYM |10] . There exist instanton solutions 
interpolating some of these vacua. Explicit instanton solutions are known. The conditions 
for allowed instanton solutions are determined by the initial and final vacua and the gauge 
group. In the representation of Young tableaux, the condition becomes moving a number 
of boxes from left to right only. The instanton actions are the difference of superpotentials 
of initial and final vacua. We compute the superpotential of each vacuum on the weak 
coupling gauge theory side, and they are traces of some operators and their final answers 
depend on the couphng and the numbers like the \ N^^'^ . 

We then study the description of these instantons in gravity duals. In the case where 
initial and final vacua are very close to each other, the instanton can be approximated as 
an Euclidean brane wrapping certain cycle in the gravity dual background. Each gravity 
background dual to the vacuum is given by the configuration of an electric charge system 
where the charges distribute on disks under an external quadratic potential. The action 
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of such Euclidean D2-brane wrapping a three-cycle shows that it is the energy difference 
between initial and final configurations of the electrostatic system. We thereby are able 
to know the instanton actions corresponding to instantons interpolating quite different 
backgrounds since we know the superpotentials. We define the superpotential of each 
vacuum on gravity side by the energy of the charge system or eigenvalue system. When 
the eigenvalues or charges do not interact with each other, they do not form disks and this 
agrees with the superpotential in weak coupling gauge theory analysis. When considering 
quantum effects of this eigenvalue system, the charges repel and form extended disks. The 
non-perturbative superpotential could be obtained by solving some integral equations. 
We also analyzed the case of a Euclidean NS5 brane wrapping a six-cycle, and it shows 
strong coupling effect of the superpotential which is not seen on the weak coupling gauge 
theory side. 

Though the instanton solutions exist, they do not contribute to the correction of the 
energies of these vacua. These instantons preserve eight supercharges and break half of 
the supersymmetries of the background. There are fermionic zero modes including the 
eight Goldstone fermions from broken supersymmetries and integrations of the fermionic 
collective coordinates in the path integral make the instanton tunneling amplitude vanish. 

We see that the geometry and the disk configuration on the gravity side is an emergent 
phenomenon from the eigenvalues of the matrices in the boundary gauge theories. On 
the gravity side, the S"^ emerges due to commutation relations of 5*0(3) scalars. The 
1] direction emerges as the eigenvalues of an adjoint scalar $. Further we study the 
emergence of and p direction qualitatively. The ground state geometries we have is 
a subsector of the 1/8 BPS geometries, associated with three i?-charges. The sector of 
the 1/8 BPS states in these theories can be described by three dimensional harmonic 
oscillators associated with the three complex scalars in these theories. They have a six 
dimensional phase space. The emerges because it is the enhanced symmetry of the 
ground state droplet in the phase space. The charges are extended so there is another 
radial direction in the phase space, which give rise to the p direction on the gravity side. 

We characterize the feature of the emergence of the electrostatic picture qualitatively. 
The forces of the charges are from the quantum repulsions of the eigenvalues. However, 
we are not able to derive the force quantitatively, which is a Green's function in three 
dimensions for the Laplace operator. In the case of the theory of U (N) A/" = 4 SYM on 
R X S^/Zj:, with k/N small, the force becomes a two dimensional Green's function, and 
this is almost the same as the TV = 4 SYM on i? x 5*^ previously known. Perhaps a more 
concrete approach to derive the force from the gauge theory side is to use a superfield 
formulation of these theories, and study the loop equations of the chiral ring operators 
in these theories, similar to the method in [66]. The loop equations might give us some 
equations that are the gauge theory results of the forces between eigenvalues. 

If we can derive the force, by the definition of the non-perturbative superpotential in 
terms of the energy of the eigenvalue system, we will be able to know the non-perturbative 
superpotentials from gauge theory side. They will also reveal the quantum nature of the 
NS5 branes in these theories. 
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The theories we studied are similar to the 3+1 dimensional = 1* SYM |25j.|26j. 
They all have many vacua and mass gaps. Some of the vacua discribe NS5 branes. Most 
of the dual geometries are smooth and similar embeddings of branes are possible if we 
know the full solutions dual to A/" = 1* SYM. The features of the gravity duals of these 
theories in this paper might help us gain some information that is similar to those of the 
J\f = 1* SYM. The NS5 branes in the M = 1* SYM wrap R^'^ x S^, while in our cases 
they wrap Rx S^, under the limit when the S"^ and become infinitely large, these two 
theories can be related via T dualities three times. This limit corresponds to expanding 
around the tip of a pair of disks in the plane-wave matrix model. There is a curious 
similarity with the eigenvalue distributions in this limit (cf. (4.2) of [H]) for the trivial 
vacuum of A/" = 1* SYM. It would be interesting to understand better the NS5 brane 
vacua and the little string theory from these theories. 

Similar to the two theories studied in this paper, the vacua of the A/" = 4 SYM on 
RxS'^Zk [10J,[27J and the IIA NS5 brane theory on i? x S*^ [lU] , [IS] are also characterized 
by charged disk configurations. It would be interesting to understand the instanton 
solutions in these theories. The disk configurations in these four theories can be related 
to each other by various limits. 

The plane-wave matrix model is also similar to the tiny graviton matrix theory in 
[65] . The model in [65] has solutions corresponding to fuzzy three-spheres and different 
vacua correspond to different representations of Spin{A). It is possible that there are 
also similar instanton solutions in [65] . This may help in understanding the instantons 
interpolating between different giant D3-branes |35j.|36]. 
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